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Abstract 

The first underground data run of the ZEPLIN-II experiment has set a limit on the 
nuclear recoil rate in the two-phase xenon detector for direct dark matter searches. In 
this paper the results from this run are converted into the limits on spin-dependent 
WIMP-proton and WIMP-neutron cross-sections. The minimum of the curve for 
WIMP-neutron cross-section corresponds to 7 x 10~ 2 pb at a WIMP mass of around 
65 GeV. 
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1 Introduction 

Weakly Interacting Massive Particles (WIMPs) remain the most plausible candidate for 
dark matter, responsible for about 80% of the total matter and for about 23% of the total 
energy contents of the Universe. WIMPs are expected to interact with ordinary matter via 
spin-independent (sometimes also called scalar or coherent) and spin-dependent (axial) 
interactions. Spin-independent interactions of WIMPs should largely dominate in high- 
A targets due to the A 2 coherence enhancement factor (here A is the atomic weight of 
the material used as a target). However, the relative probability of the spin-independent 
and spin-dependent interactions depends also on the particle content of WIMPs, i.e. on 
parameters of a particular supersymmetric model. In some models, WIMPs interact pre- 
dominantly through the spin-dependent interactions. This stimulates the search for spin- 
dependent interactions of WIMPs in addition to the spin-independent case. 

The interest in spin-dependent WIMP-nucleon interactions is enhanced by the positive 
WIMP signal claimed by the DAMA Collaboration. Although the first DAMA publica- 
tions favoured spin-independent signal due to the very soft spectrum of events observed, 
further studies allowed the interpretation of the results in terms of the combination of 
both spin- independent and spin-dependent interactions [T]. 

The spin-dependent WIMP-nucleus cross-section depends on the spin factor of the nucleus 
that is primarily determined by the number of protons and neutrons in the nucleus, namely 
whether it is odd or even. For odd-proton nuclei the spin-dependent WIMP-nucleus cross- 
section is mainly due to the WIMP-proton interactions, whereas for odd-neutron nuclei 
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it is dominated by WIMP-neutron scattering. For even-even nuclei the spin-dependent 
cross-section is highly suppressed. 

Among the most sensitive targets for spin- dependent WIMP searches are Na, I, Cs, F 
(all odd-proton) and Xe, Ge (having odd-neutron isotopes). The best limits from direct 
detection experiments so far have been set by the NAIAD (WIMP-proton) [2] and CDMS 
[3], ZEPLIN-I [4] and EDELWEISS [5] (WIMP-neutron) experiments. 

In this paper, the recent results of the ZEPLIN-II experiment are analysed in terms of the 
spin-dependent limits. A brief description of the detector and experimental data is given 
in Section 2. The method of calculating spin-dependent limits is presented in Section 3 
together with the results. Conclusions are given in Section 4. 



2 ZEPLIN-II experiment 

ZEPLIN-II [6|f7] is a two-phase (liquid/gas) xenon detector searching for elastic scat- 
tering of WIMPs off xenon nuclei. ZEPLIN-II is operated at the Boulby Underground 
Laboratory in the UK at a depth of 2805 m w.e. underground with a muon flux of 
(4.09±0.15) xlO -8 muons/cm 2 /s [8]. The detector, data acquisition system, analysis pro- 
cedure and experimental data are described in detail in Refs. [9)10)11] , We present here 
only a short summary of the detector performance and some other features important for 
further analysis. 

ZEPLIN-II consists of a vacuum cryostat containing about 31 kg of liquid xenon. The tar- 
get volume is viewed from above by seven 13 cm diameter photomultiplier tubes (PMTs). 
The inner surface of the xenon vessel is covered by light-reflecting PTFE for better col- 
lection of scintillation photons. 

The detector records scintillation light and ionisation electrons from charged particles. 
Scintillation light is detected by the PMTs promptly after excitation of the active medium. 
Ionisation electrons are drifted through the liquid towards the surface by means of applied 
electric field, extracted from the liquid into the gas where they are accelerated producing 
a secondary scintillation signal as electroluminescence. The depth of liquid xenon (14 cm) 
corresponds to a maximum drift time of 73 /as for electrons at the drift field of 1 kV/cm. 

To protect the target volume from radioactivity in rock (gammas and neutrons), the 
ZEPLIN-II detector is surrounded by hydrocarbon material and high-purity lead. Part of 
the hydrocarbon shielding is instrumented liquid scintillator that also acts as an anticoin- 
cidence system (active veto) preventing the signals, detected simultaneously in the target 
and in the veto, being interpreted as WIMP interactions. 

The trigger is provided by five-fold coincidences between different PMTs at a single pho- 
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toelectron level. Either a primary, SI, or a secondary, S2, signal can trigger the data 
acquisition. The signals from all seven PMTs are recorded with 2 ns sampling time during 
±100 /is around the trigger pulse, covering all possible arrival times for both primary and 
secondary signals. In an off-line data analysis three-fold coincidences between different 
PMTs at a single photoelectron level were used to identify and parameterise the primary 
signal SI. 

A number of parameters have been measured for each waveform. They are listed and 
discussed in Refs. [9fT0] . the most important being the total areas of the primary and 
secondary pulses (proportional to the number of photoelectrons) , SI and S2, respectively, 
the time delay between them and the width of the pulses (that determines whether the 
pulse is the primary or secondary signal). 

Position sensitivity of the experiment in the vertical direction is achieved by considering 
the time delay of the secondary pulse which is proportional to the drift time of the 
electrons and thus determines the point of the interaction along the drift field direction. 
This allows us to exclude events originating on or close to the grid wires that provide the 
electric field and are contaminated mainly with radon progeny. In the horizontal plane 
the event position is reconstructed using the relative pulse areas from secondary signals 
on different PMTs. This method gives bigger uncertainty compared to the drift time due 
to the large PMT sizes and small photon statistics at low energies. 

A daily energy calibration of the detector using a 57 Co gamma-ray source allowed moni- 
toring the stability of detector operation. The WIMP / gamma discrimination performance 
of ZEPLIN-II has been tested by calibrating the detector using high-energy gamma-ray 
( 60 Co) and neutron (AmBe) sources. High-energy gamma-rays produce the main electron 
recoil background, whereas fast neutrons scatter elastically off nuclei producing nuclear 
recoils in the same way as expected from WIMP scatters. Using the S2/S1 versus SI 
plot from neutron calibration run (Figure 8 in Ref. [§]) the nuclear recoil acceptance box 
has been defined as retaining 50% of nuclear recoil events at any given energy chosen for 
analysis. This acceptance box was expected to have a small number of electron recoils due 
to the tail of S2/S1 distribution observed in the gamma calibration run. 

The first data run of the ZEPLIN-II detector had a live time of 31.2 days after time 
periods with unstable operation conditions were removed from the analysis. A number 
of software cuts have been applied to the parameterised pulses allowing the selection of 
single interactions within the fiducial volume of the detector to be made. These cuts are 
described in Refs. [9]. For each cut an energy dependent efficiency has been evaluated 
either from the data or from a combination of data and simulations. 

An important cut that reduces significantly the fiducial volume of xenon, is the radial 
cut. This rejects the events that are reconstructed as being close to the PTFE walls. In 
reality, due to imperfect position reconstruction of the interaction points in the horizontal 
plane, a long tail of events assumed to be originated at the walls is reconstructed towards 
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the centre of the detector. The majority of these events is believed to be caused by the 
alpha decay of radon progeny accumulated on the charged surfaces. To remove most of 
the 'wall' events the radial cut has been applied reducing the fiducial mass of xenon down 
to 7.2 kg (see Ref. [9] for more details). 



Two xenon isotopes occurring naturally are sensitive to WIMP-nucleus spin-dependent 
interactions: 129 Xe and 131 Xe. Other stable xenon isotopes are even-even (with even num- 
bers of proton and neutrons) with very small coupling to WIMP spin. 129 Xe and 131 Xe are 
both odd-neutron isotopes and hence are sensitive mainly to WIMP-neutron interactions. 
Their relative abundances are 26.4% and 21.2%, respectively. The effective exposure of 
these target isotopes to WIMPs is calculated as Xj = M x t x d x A x e, where Xi is the 
exposure of the ith isotope, M = 7.2 kg is the fiducial mass of the target, t = 31.2 days 
is the live time of the run, is the relative abundance of the isotope, A = 0.5 is the 
fraction of nuclear recoils in the acceptance window on S2/S1 vs SI plane chosen for the 
analysis and e is the overall energy dependent efficiency of other cuts described above and 
in Ref. [9]. If we neglect the energy dependent efficiency e which increases from about 
35% at 5 keV to about 75% at 10 keV, the exposures for different xenon isotopes are 
29.7 kgxdays for 129 Xe and 23.8 kgxdays for 131 Xe. Cut efficiency e further reduces the 
exposure making it also energy dependent. 



The energy range of 5-20 keV (electron equivalent) has been chosen for the data analysis. 
Below 5 keV the trigger efficiency is rather small (less than 40%). Above 20 keV the 
sensitivity of xenon target to WIMP interactions decreases significantly because of the 
rapidly falling form-factor. To avoid any bias, the selection of cuts and energy range for 
the final analysis was based on the results from calibration runs and from 'unblinded' 10% 
of data prior to 'opening the box' with the remaining 90% of data. 



The analysis of data has revealed 29 events in the nuclear recoil acceptance box, the 
expected background rate due to electron recoils and 'wall' events being 28.6 ± 4.3 (see 
[9] for full description of the procedure to evaluate the expected background). Based on 
the previously published simulations [T3|14|15p6] we expect to have less than 1 nuclear 
recoil from neutron background in the detector for the aforementioned exposure. 



Applying the procedure described by Feldman and Cousins [T7] the 90% CL upper limit 
on the number of nuclear recoils has been set as 10.4 using the ROOT software [32] and 
then converted into an upper limit on the WIMP-nucleon spin-independent cross-section 
[9]. The limit on the nuclear recoil rate can also be used to set a limit on spin- dependent 
interactions. This procedure is described below. 
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3 Limits on spin-dependent cross-sections 



In the spin-independent case Majorana WIMP coupling to protons and neutrons is very 
similar, the coherence enhancement factor is proportional to A 2 and the cross-section on 
the nucleus does not depend on particular WIMP model parameters. For spin-dependent 
interactions the coupling to protons and neutrons is very much different and depends 
strongly on the WIMP model parameters. In the derivation of the spin-dependent limits 
we follow the procedure described in Ref. [15) . 

The WIMP-nucleus cross-section, a a, can be written as: 

a A = — G F fi A {a p < b p > +a n < S n >) — — , (!) 



where Gf is the Fermi weak coupling constant, \ia is the WIMP-nucleus reduced mass, a V)Tl 
are the effective WIMP-proton (neutron) couplings, < S Pjn > are the expectation values 
of the proton and neutron spins in the nucleus (or spin-factors) and J is the nuclear spin. 
For a proton or a neutron as a target, Eq. (CQ) is transformed into the cross-section for 
WIMP-proton (neutron) interactions with the proton (neutron) spins < S Pjn >= 1/2 and 
J = 1/2. Eq. ([1]) does not correspond to the total WIMP-nucleus cross-section, but to that 
at zero momentum transfer. As in the case for spin-independent interactions, the cross- 
section for zero momentum transfer is usually presented for comparison with other results 
and model predictions. To evaluate the 'real' interaction cross-section and to compare it 
with the experimental data this has to be multiplied by the nuclear form-factor F 2 (q) 
which is a function of the momentum transfer, q [20J. 

As in Ref. [19] we assume that the total WIMP-nucleus cross-section at zero momentum 
transfer is dominated by either WIMP-proton or WIMP-neutron interactions only, setting 
the 2nd component equal to 0. In this case the WIMP model dependent parameters are 
cancelled in the equation for the ratio of the cross-sections: 

^7 ~ 174 <s P , n > 2 J + i' 1 ' 



which retains only the nuclear physics parameters. As a a is measured in an experiment 
or a limit on a a is set from the experimental data (assuming a particular form-factor) 
and nuclear physics parameters < S p>n > can be calculated independently of the WIMP 
model, the WIMP-proton (neutron) cross-section can be evaluated in a (almost) model- 
independent way using Eq. ([21 . The only WIMP model dependence that remains is hidden 
in the form factor that appears to be different for different WIMP models. This depen- 
dence, however, is not significant for most isotopes and is comparable to the uncertainty 
in the nuclear physics model used. 
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For xenon as a target, additional complication comes from the existence of two isotopes 
with odd-neutron nuclei, 129 Xe and 131 Xe. In this case the combined limit on the cross- 
section is calculated as [20J: 



/ i i 

P' n A, a p,n 



where a^ n are the WIMP-proton (neutron) limits set for cross-sections on individual 
isotopes taking into account their fraction by weight, and a p , n is the combined limit on 
the WIMP-proton (neutron) cross- sect ion. 

In further calculations we used the spherical isothermal dark matter halo model with the 
following parameters: pd m = 0.3 GeV cm -3 , v = 220 km/s, v esc = 600 km/s and VEarth = 
232 km/s. The form-factors were computed for the two xenon isotopes using the nuclear 
shell model calculations [2TJ with 'Bonn A' nucleon-nucleon potential. This is based on a 
comprehensive meson-exchange model for the nucleon-nucleon interaction in field theory 
developed by the Bonn group [22]. This model offers the most consistent approach to the 
nuclear many-body problem at low energies relevant to WIMP interactions. It includes all 
important diagrams with a total exchanged mass up to about the cutoff mass (~ 1 GeV). 
The various meson-exchange contributions in this mass range are introduced step by step 
proceeding from lowest-order to higher-order processes and from long range to short range. 
The model predictions agree well with the deuteron data and parameters derived from 
nucleon-nucleon scattering experiments [22]. 

The total form-factor can be written in the form: 

F 2 () = , A) 

W 5(0)' K ' 

where 

S(q) = alS 00 (q) + a\S n {q) + a^S^q). (5) 



Here a = a p +a n , a\ = a p —a n and Sqo, Su, Soi are the isoscalar, isovector and interference 
contributions to the spin structure function S(q), respectively. These contributions are 
determined by a nuclear model (independent of the WIMP type), with coefficients ao and 
ai related to the WIMP-proton and WIMP-neutron coupling constants a p n (and hence 
WIMP- type dependent). The coupling coefficients depend also on the assumption about 
the quark spin fractions, i.e. the fractional contributions of different quark species to the 
nucleon spin. We used the WIMP-proton and WIMP-neutron coupling constants a p and 
a n with the quark spin fractions from Ref. [23]. Although the form-factor is normalised 
to the value at zero momentum transfer (Eq. (J3])) that also includes the WIMP coupling 
coefficients, they do not cancel out completely, leaving some dependence on the WIMP 
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particle model. The dependence on the WIMP model parameters for xenon isotopes is 
not as weak as, for instance, for iodine and low-A isotopes. The form-factor for higgsino 
interactions, however, is the smallest among WIMP particle models, so the limit based 
on the present calculations with higgsino form-factor is conservative. 

The spin-factors < S PtTl > for the two odd-neutron Xe isotopes have been calculated in 
Ref. [21J and also given in Ref. [19]. In the present analysis we used the values reported 
for the same nuclear model with Bonn A potential as for the form-factors. 

Figure [1] shows the 90% CL limits on the WIMP-proton (a) and WIMP-neutron (b) 
spin-dependent cross-sections calculated using the ZEPLIN-II data and the procedure 
described above. Results from some other experiments are also shown. The ZEPLIN-II 
limits are also given in numerical form in Table [1] to allow more accurate comparison 
with other experiments. The minimum of the curve for the WIMP-neutron cross-section 
corresponds to 7 x 10~ 2 pb. The limits are dominated by the contribution from the 129 Xe 
isotope. The ZEPLIN-II limits on WIMP-neutron cross-section are comparable to the 
currently best result obtained by the CDMS experiment [3|. Although the CDMS spin- 
independent limits are better than those from ZEPLIN-II, the spin-dependent WIMP- 
neutron cross-section limits presented here are very similar to those set by CDMS due to 
the higher fraction of odd-nucleon isotopes in xenon. 

The uncertainties in the nuclear spin and form-factors used for the evaluation of limit 
are not negligible. Apart from the remaining dependence on the particle model, there 
is an uncertainty related to the nuclear model. Ressell and Dean [21] found a factor of 
2 difference between their calculations of the spin structure functions for 131 Xe (higher 
values at zero momentum transfer) and earlier calculations by Engel [27J. The model 
used by Ressell and Dean [21 J with Bonn A potential gives smaller (again by a factor 
of 2 approximately) spin structure functions than simple 'single particle' model. Slightly 
lower spin expectation value S n (by about 20%) for 129 Xe was obtained with Nijmegen 
II potential [21J which results in a 20% higher limit on the cross-section. This makes 
the systematic uncertainty of the cross-section limit as large as a factor of 2 due to the 
nuclear model calculations. We stress, however, that the model used in the present analysis 
is based on the most recent and accurate calculations of nuclear spin and form-factors 

m- 

The predictions of different nuclear models for different nuclei are not strongly correlated 
in the sense that two models can predict, for instance, similar spin expectation values 
for one nucleus but very much different spins for another nucleus. Thus, possible future 
detection of WIMP interactions with several different target nuclei should help in reducing 
the uncertainties in nuclear physics models and, hence, in improving accuracy of the 
WIMP parameters' estimates. 

The limits shown in Figure [T] were obtained for a pure higgsino as a WIMP. Assuming 
pure photino or bino as a WIMP gives 35-40% better limit on the cross-section, whereas 
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the limit for pure zino would be a factor of 3 lower. As precise composition of WIMPs is 
not known, the assumption of a pure higgsino leads to the most conservative limit. 

To assess the significance of the spin-dependent limits less than 0.1 pb, it should be 
noted that a marginal dark matter candidate, the heavy Majorana neutrino, could have 
a spin-dependent cross-section of about 0.01 pb. However, a re-evaluation of likely super- 
symmetric candidates by Ellis [28] suggests spin-dependent cross-sections no greater than 
10 1 pb. 

Figure [2] shows constraints on the WIMP-proton and WIMP-neutron coupling coefficients 
a p and a n (for a WIMP mass of 50 GeV) as derived from the ZEPLIN-II data in comparison 
with other experiments. Here again we used the formalism described in Ref. [19] that 
allows conversion of the cross-section limits into the allowed regions on the a p — a n plane 
using the equation: 



where the small mass difference between the proton and the neutron is ignored. 

Strictly speaking the procedure is mathematically correct only if the form-factor used to 
evaluate the limits on cross-sections a p l n is independent of the WIMP model. In most 
cases (including ZEPLIN-II limit derivation in this paper) the form-factor depends on 
the WIMP model, i.e. on the WIMP-proton and WIMP-neutron coupling coefficients 
a p>n . In practice this means that the limits a p ^ n depend on a Ptn . We used the form-factor 
for pure higgsino and in the derivation of the allowed regions on the a p — a n plane we 
ignored the dependence of the form-factor, defined as F 2 (q) = S(q)/S(0), on the coupling 
coefficients. This remark concerns also other experiments [2,5,24,25] for which similar 
procedure was used. We repeated the derivation of the allowed regions on the a p — a n 
plane from the cross-section limits from these experiments and found them to be in good 
agreement with the original publications [5,2 4"|[2"5] . The allowed region from the CDMS 
experiment was copied from Ref. [3J. Our computation of the allowed region for CDMS 
gave slightly different result which implied that the CDMS Collaboration used another 
procedure for constraining coupling coefficients in Ref. [3j. The allowed region from the 
DAMA/Nal experiment [1] was calculated using the allowed regions for the cross-sections 
from Ref. |26j . 

One of the alternative methods of setting limits on the coupling constants was suggested 
in Ref. [26]. Despite its complexity, this method has an advantage of providing model 
independent constraints on the coupling coefficients by using the spin structure function 
S(q) directly in the process of the data analysis, i.e. the coupling coefficients a p>n being free 
parameters in the fit to the data. Note, however, that for accurate comparison between 
different experiments, all data have to be analysed using the same method. 
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Figure [2] shows that for 50 Ge V WIMP mass the DAMA allowed region (filled area on the 
figure) is excluded by the combination of other experiments. Savage et al. [26], however, 
found that the interpretation of the DAMA positive signal in terms of spin-dependent 
interactions is still compatible with other experiments at small WIMP masses (5-13 GeV). 
NAIAD (sensitive mainly to a p ) and CDMS/ZEPLIN-II (sensitive mainly to a n ) provided 
so far the most stringent constraints on the WIMP-nucleon spin-dependent interactions. 



4 Conclusions 

The upper limits on the WIMP-proton and WIMP-neutron spin-dependent cross-sections 
have been set using the ZEPLIN-II data. The minimum of the curve for WIMP-neutron 
cross-section corresponds to 7 x 10 -2 pb. The limits on WIMP-neutron cross-section are 
comparable to the currently best result obtained by the CDMS experiment. 
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Table 1 

90% CL limits on spin-dependent WIMP-neutron, a n , and WIMP-proton, a p , cross-sections as 
functions of WIMP mass, M\y, for 129 Xe and 131 Xe isotopes together with the combined limits 
from the ZEPLIN-II experiment. 
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Fig. 1. 90% CL upper limits on the WIMP-proton (a) and WIMP-neutron (b) spin-dependent 
cross-sections derived from ZEPLIN-II data. The limits from the two isotopes ( 129 Xe - dashed 
curves, 131 Xe - dotted curves) and the combined limit (solid curves) are presented. The limits 
from some other experiments are also shown: NAIAD [2] (•), ZEPLIN-I [4J (■), CDMS [3] (o), 
EDELWEISS [5| (□) and PICASSO [Mj (A) (the latter result coincides with that of the SIMPLE 
experiment [25]). The interpretation of the positive annual modulation signal observed by the 
DAMA experiment [I], in terms of the constraints on spin-dependent cross-sections reported by 
Savage et al. [26], is shown by the filled area. 
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Fig. 2. Constraints on the WIMP-proton and WIMP-neutron coupling coefficients a p and a n 
(for a WIMP mass of 50 GeV) as derived from 90% CL upper limits on the WIMP-proton and 
WIMP-neutron cross-sections from different experiments. The ZEPLIN-II results are shown by 
solid curves (almost vertical parallel lines). The region between the two lines is allowed by the 
ZEPLIN-II data. Other results are: CDMS (dashed curves) f3j, EDELWEISS (dotted curves) 
[5], NAIAD (dashed-dotted curves, nearly horizontal) [2j and PICASSO/SIMPLE (solid curves, 
nearly horizontal) |24|25j . The allowed regions from the DAMA/Nal experiment [T] is shown by 
the filled area as calculated from the cross-section allowed regions reported in Ref. [26] . 
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